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Abstract 

Despite suggestions to the contrary, we show in this paper that the usual dispersive form of the 
electromagnetic energy must be used to derive the Lifshitz force between parallel dielectric media. 
This conclusion follows from the general form of the quantum vacuum energy, which is the basis 
of the multiple-scattering formalism. As an illustration, we explicitly derive the Lifshitz formula 
for the interaction between parallel dielectric semispaces, including dispersion, starting from the 
expression for the total energy of the system. The issues of constancy of the energy between parallel 
plates and of the observability of electrostrictive forces are briefly addressed. 
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I. INTRODUCTION 



Recent years have yielded considerable progress in understanding quantum vacuum or 
Casimir energies, both theoretically and experimentally. For a very recent review see Ref. 
However, there are controversial aspects, both having to do with the concept of zero-point 
energy applied to a single system, or to the universe as a whole and with including 
thermal corrections, and their observability in experiment The latter question refers to 
how the electric properties of materials depend on (imaginary) frequencies, that is, upon 
dispersion. 

In this paper, we address the latter issue. In a recent paper we had proposed, following 
a suggestion of Lifshitz j4|, that the usual dispersive term in the electromagnetic energy for 
a given frequency jst , 



-E\y) + H\y) 



(we ignore the magnetic susceptibility, that is, we set = 1), should not be included. 
However, the usual derivations of the Lifshitz interaction between dielectric slabs are not 
based on the total energy. For example, in Ref. {g] the Lifshitz formula is derived from 
the pressure, or equivalently the spatial components of the stress tensor, and also from the 
variational principle enunciated in Ref. |7j. It is also easy to obtain this sa^e^result using the 



3. Equivalently, 



recently repopularized multiple-scattering approach to Casimir energies 
the multiple- reflection expansion yields the Lifshitz formula immediately 10 1. 

In this note, we derive the Lifshitz energy directly from Eq. (11. ip . We first see, in 
Sec. [ni how dispersion is incorporated in a general formulation. This demonstrates that the 
dispersive form of the energy is required. Then, after giving the form of the Green's dyadic 
in Sec. IIIIl in Sec. |IV]we will explicitly derive the Lifshitz formula from Eq. (II. ip . and will 
see manifestly that the dispersive term provides the Jacobian of the required transformation 
of coordinates necessary to obtain the necessary logdet form. In the Conclusions, we also 
bring up the related possibility of measuring electrostrictive effects in liquids. The Appendix 
points out that the well-known constancy of the energy density between parallel perfectly 
conducting plates does not hold for dielectric plates (that is, if regions 1 and 2, defined in 
Sec. mil are constituted of dielectric material), or even if a dispersive medium exists between 
metallic plates. 
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II. GENERAL FORMULATION 



Let us start from Eq. (11 .ip . and consider the quantum vacuum energy associated with 
electromagnetic field fiuctuations: 



e = - I (dv) 



du 
2^ 



duo 



(2.1) 



The expectation values appearing here are given by the electromagnetic Green's dyadic, 



(E(r)E(r')) 
(H(r)H(r')) 



-r(r,r'), 
1 1 



V X r(r,r') X V', 



(2.2a) 
(2.2b) 



and so inserting these into the energy expression (12. ip . integrating by parts, and using the 
differential equation satisfied by the Green's dyadic. 



V X V X r + Lo^eT 



we obtain the expression for the energy 



dco 
2^ 



de 

2e trF + w— trF 

du 



(2.3) 



(2.4) 



We can obtain this same result starting from the standard trace-log formula: 
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2 / 27r du 2 2-n du 



2 J 2n 
where from Eq. (12. 3p 
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-V X V X -e, 



u^ 



(2.5) 



(2.6) 



and where Tr includes the trace over spatial coordinates. The final form in Eq. (12. 5 p is 
exactly the result (12. 4p derived from the expectation value of the classical electromagnetic 
energy (II. ip . 

To conclusively demonstrate that the dispersive term must be included, we derive the 
variational principle used to obtain the Lifshitz formula in Ref. 3]- This depends upon the 
variational statement 

5F = -F(5F~^F = FfeF. (2.7) 



Also, using the differential equation for the Green's dyadic we find 



dT 



r— — r = -TeT + r— r + -r. 

dcu CO du) u) 



Therefore, the e-variation of Eq. (12 ■4p yields 

2 J ZTT \ du du J 

= -U ^u^^Tr 6eT, 
2 J 271 du 

which, upon integration by parts, yields the variational principle used in Refs. 

du 



(2.8) 



(2.9) 
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See also Ref, 
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III. GREEN'S DYADIC FOR PARALLEL SLABS 



In this and the following section we supply an explicit derivation of the Casimir-Lifshitz 
interaction between parallel dielectric slabs (of infinite thickness). Specifically, consider a 
dielectric function in the following form 



ei, 2; < 0, 
e(r] = < 63, < z < a, 
€2, a < z. 

Then the Green's dyadic can be written as a transverse Fourier transform. 



rfr,r') 



{dk± 
T2^ 



^^±i'^-'^')±^(z,z';k^,Uj, 



where the reduced Green's dyadic has the form 



B9, 



e dz e' dz'" ee' dz" 



u^g"" 

ik d „E n 



(3.1) 



(3.2) 



(3.3) 



Here we have dropped 5-function terms, we have denoted e = e{z), e' = e{z'), and we have 
chosen the coordinate system so that k_L = kct. Here the TE {H) and TM {E) (relative to 



4 



the z axis) Green's functions satisfy the differential equations 

92 



d 1 d 



dz edz e 



5{z-z'), 
6{z-z'). 



(3.4a) 
(3.4b) 



We will solve these equations in each of the three regions given in Eq. (13. ip . subject to 
boundary conditions between the regions that and dzQ^ are continuous, and that and 
{l/e)dzg^ are continuous. These boundary conditions reflect the underlying requirement 
that the transverse parts of E and H are continuous, while the normal component of D = eE 
is continuous (there are no surface charges or currents). It is a straightforward calculation 
to find the Green's functions in each region. We display the results for the only situation we 
need in the following, when z and z' are both in the same regions. Below the first interface, 



z, z' < 0, 



ri 



1 

Ki + 



+ 



4kiK3 1 



(3.5a) 
(3.5b) 



where = k"^ + ("^EaiiC), a = 1,2, 3, and we have made a Euclidean rotation, u = i(. Here 
we have introduced the abbreviation 



d 



(3.6) 



Similarly, above the second interface, z, z' > a. 



g''{z,z' 



r2 



1 

K2 - K,3 



^-K2\z~z'\ _|_ ^^g-K2(2+z'-2a) 



f^2 + /t3 

In the intermediate region, a > z, z' > 



+ 



4K2/t3 1 



(3.7a) 
(3.7b) 



9''{z,z') 



1 
2^ 



e "^l"' "^'1 + -coshK3(z - z' 
d 



1^3 + ^2 l^-K3(2+2'-2a) 
— K2d 



(3.8) 



^3 — /?! d 

The transverse magnetic Green's function g^ is obtained from the above by replacing 
f^a f^a/^a Gxccpt in the exponents. 
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The fact that g{z,z) depends on z imphes, in general, that the mean-squared electric 
and magnetic fields also depend on position, as does the energy density. This seems to 
contradict the fact that for parallel conducting plates the energy density is constant in each 
region js], We shall show, in fact, in the Appendix that the electromagnetic energy 
between perfectly conducting plates is indeed constant, provided the intervening medium is 
nondispersive. 



IV. LIFSHITZ ENERGY 

The Casimir-Lifshitz energyper unit area for the situation of parallel slabs described by 
the dielectric function (13. ip is |6| 



1 POO POO 

^=T-2 dkk[hi{l- rTET^^e-^'""') + In (l - rTMrUe"'''^")] 



(4.1) 



where 



rr,^^, = (4.2) 

and the TM reflection coefficients are obtained by replacing Ka ^ i^'a = i^a/^a- In this 
section, we rederive this result from Eq. (12. 4p . 



A. TE contribution to the energy 

The TE part of the energy can be written as 

1 poo roo POO 1 

^™ = ^y^ dCj^ kdkj dzC^{-K')g^{z,z). (4.3) 

In each region, the dispersive term is necessary to change variables from ^ to k^. In the first 
region, omitting infinite terms which contain no reference to the separation a between the 
regions, we find rather immediately (we assume, as usual, C^ea(C) — as C — > 0) 

1 /-OO POO ^ 

'^™'' = -7r^/ f'dk InATE, (4.4) 

where the partial derivative means that k,2 and K3 are not altered, and 

Ate = 1 - r-TEr^Ee^'"'"- (4-5) 
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Similarly, in region 2, 



= -TTT / kdk dK2^—\nATE. (4.6) 



The intermediate region involves a slightly more involved calculation, but the result has the 
same form (after we omit a constant term in the force): 



where now the derivative acts also on the exponent in Ate- In this way we obtain exactly 
the expected TE contribution: 

^TE L r kdk r dn ^ (— + —— + In A 

2vr2 7o Jk ^2\dKi dKsdK2 dnsdnij ™ 

POO 

dkk rfClnATE, (4.8) 
Jo 



47]-2 

which is just the first term in Eq. fl4.ip . 
B. TM contribution to energy 

The TM contribution requires a somewhat more elaborate calculation. The TM contri- 
bution to the trace of the Green's dyadic is 

trg^ = --^^^g^ + ^g^. (4.9) 

e oz e' oz' ee' 

This differential structure has different forms depending on whether it acts on the pure 
exponential terms in g^ , or on the hyperbolic cosine in Eq. (13. Sp . namely, in the first case, 

I d 1 d k"^ I , .0 ^ .X 

and in the second case, 

I d I d k"" C 

— H ^ - — . (4.11) 

e oz e' oz' ee' 63 

Except for that last exceptional case, combining this trace term with the dispersive term in 
the energy gives 
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Thus, in region 1, the contribution to the TM energy is 

^TM,l 



(4.13) 



where Atm differs from Ate by replacing Ha by k'^ = Ha/^a except in the exponents, that 
is, tte — ^ '^TM- Similarly, in region 2, 



TM,2 



(4,14) 



In region 3, however, we have to take into account the special case (14. lip . It is most 
convenient then to regard ^3 and k'^ as independent, in which case we can write 



bTM,3 



1 



47r2 



00 fOO 

k dk 



dC 



'0 ^0 
Thus, the total TM contribution is 



^3 - - 



TM 



47r2 



dkk\ \ dC 

{Jo 
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2 dC J dn'^ 2 dC dn^ 



In A 



TM- 



(4.15) 



Sk'^ L/n.2 L/n.3 

C f d dn'^ d dK'2 d dn'^ d 



+ K 



d 



dK'i — 



+ 



In Atm 

In Atm 



(4.16) 



r ' 2 \ (9^3 dn^ dK^ dK'2 dK3 dK[ 

The first term here is actually zero, because the differential operator annihilates Atm? since 



J 9 



and so we obtain exactly the Lifshitz result 

1 



^l^^TM = -«3^^TM 



(4.17) 



TM 



47r2 



kdk (i^lnATM- 



(4.18) 



V. CONCLUSIONS 



Ordinarily one calculates the Casimir-Lifshitz free energy directly from the pressure, or 
from an equivalent variational approach. Therefore, it was not obvious how the dispersive 
term present in the energy in order to have the required balance between energy and mo- 
mentum, as in the electromagnetic energy-momentum tensor, plays a role. Earlier we had 
suggested sj] that such a term simply be omitted. However, we now see that the dispersive 
term is precisely what is needed to achieve agreement between the different formulations of 
the energy, and that the dispersive term provides the Jacobian factor necessary to derive 
the Lifshitz free energy from the expectation value of the electromagnetic energy. 
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The following point, related to the possibilities of experimental observations, ought to be 
noticed: As we have seen, a characteristic property of dispersion is that the factor d{eijj)/du 
occurs in the energy and not in the pressure or the stress. This has a bearing on the famous 
Abraham-Minkowski energy-momentum problem. As is known, an important experiment in 
this area is the Jones- Richards radiation pressure experiment showing how the effective 
pressure against a mirror immersed in a liquid varies with respect to the refractive index (cf. 

n n 

also the follow-up experiment of Jones and Leslie [ij]). The book by Jones [15[ contains 
a nice exposition of these very accurate experiments. The electrostrictive forces do not 
contribute to the radiation pressure. 

Does this imply that electrostrictive forces in a liquid are generally non-observable? Not 
quite so, although a difficulty is that at thermal equilibrium the electrostrictive forces give 
rise to elastic pressures in the liquid, acting in the opposite direction. There are ways to over- 
come this difficulty, however. One option is to proceed as in the Goetz-Zahn non-equilibrium 



experiment 
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17l |: cf. also the detailed discussion on this experiment in Ref. [18||, p. 149. 
One applies an electric field with high frequency u between two condenser plates in a liquid, 
and measures the attractive force between the plates for instance by means of a piezoelectric 
transducer. The point is that u must be so high that the elastic pressure does not have time 
to built itself up. The critical parameter here is thus the velocity of sound in the liquid. 
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Appendix A: Constancy of energy for conducting plates 

Consider the case of parallel perfect conducting plates separated by a nondispersive 
medium with dielectric constant e. The TE Green's function is [obtained by taking ki^2 oo 
in Eq. 

H_ 1 f ^_.|,_,.| , 2 cosh Kjz - z') - e-(-+-') - e-(-+-'-^") ) 
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where k'^ = k"^ + ("^e. The TE energy density is given by 



where 



, 27r J (27r)2 
Then we easily see that 



u 



TE 



dC 



-oo 271" JO 



dkkl 1 



2tt k e^'^ 



2Ka 



'2K{z—a)'\ 



The TM Green's function for perfectly conducting plates has a similar form, 
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1 f ,n 2 cosh k{z - z') + e^(^+^') + £-''(^+^'-2'^) 



2 k 



o2Ka 



(A2) 

(A3a) 
(A3b) 

(A4) 

(A5) 



and the corresponding energy density is 

1 ,„o „o, 1 



u 



TM 



27r y (27r)2 
1 



(9A' + A:V(^,^')L, 



2e 



(iC dkkl 
Jo 



271 L 271 K e'^'"' - 1 



A-2 



(A6) 



The z-dependent terms exactly cancel between Eqs. (IA4p and (IA6P and the remaining 
terms are equal, and sum to the usual Casimir energy density. 



u 



Svr^v^ Jo 



dn K 



7T 



o2K,a 



720 V^a^ ■ 



(A7) 



This cancellation, resulting in the constancy of the energy density, is rather special, 
however. It does not occur if dispersion is present, de/d( ^ 0, in which case the local energy 
density has the nonconstant form: 



u 



dC f {dk±) 1 1 
2^ 



I 2 + -—I + ('p2«^ + p-2K(2-a)^ 



) . (A8) 



Nor can the cancellation occur for dielectric media constituting regions 1 and 2, since the TE 
and TM reflection coefficients are then different. Nevertheless, we note that the z integral 
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of the spatially varying part of the energy is a constant, independent of a, and so does not 
contribute to the force on the plates. This is just as occurs for a nonconformally coupled 
massless scalar field confined between Dirichlet plates. 
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